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Abstract 



Using the reduced density matrix formalism the time dependence of 
the exciton scattering in hght-harvesting ring systems of purple bacte- 
ria is calculated. In contrast to the work of Kumble and Hochstrasser 
(J. Chem. Phys. 109 (1998) 855) static disorder (fluctuations of the 
site energies) as well as dynamic disorder (dissipation) is taken into 
account. For the description of dissipation wc use Rcdficld theory in 
exciton eigcnstatcs without secular approximation. This is shown to be 
equivalent to the Markovian limit of Capck's theory in local states. Cir- 
cular aggregates with 18 pigments are studied to model the B850 ring 
of bactcriochlorophyls within LH2 complexes. It can be demonstrated 
that the dissipation is important for the time-dependent anisotropy of 
the fluorescence. Smaller values of static disorder are sufficient to pro- 
duce the same decay rates in the anisotropy in comparison with the 
results by Kumble and Hochstrasser. 
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1 Introduction 



The physical motivation of the present paper is a growing interest in the exci- 
ton transfer in the antenna systems (ASs) of the purple bacteria photosynthetic 
units (PSUs). The light-harvesting complexes that are present in photosyn- 
thetic systems perform two major functions: harvest (absorb) the incident 
light and transport the energy in form of Frenkel excitons to the reaction cen- 
ter (RC). Many investigations, both experimental and theoretical ones, have 
been directed towards understanding of the exciton transfer in the ASs of pur- 
ple bacteria PSUs and are reviewed in and references therein]. Although 
considerable progress has been made during recent years, our knowledge about 
the mechanism of energy transfer and relaxation is still far from complete. 

In 1995 the first high-resolution three-dimensional X-ray structure of a bac- 
terial antenna complex LH2 from Rhodopseudomonas acidophila was published 
by Mc Dermott et al. @,|^. It consists (see Fig. 2 in 0) of nine identical units, 
protomers, each of which is formed by two proteins {a and (3) with bound bac- 
teriochlorophyls (BChl). The nine a subunits are packed in an inner ring to 
form a hollow cylinder of radius 1.8 nm. The 9 (3 subunits are arranged radially 
outwards with respect to the a subunits to form another ring with a radius 
of 3.4 nm. The protein serves as a scaffold for the BChls and furthermore 
specifically influences the spectroscopic properties of the BChls by supplying 
a characteristic environment for them. A ring with a radius of about 2.5 nm 
of 18 BChl molecules is sandwiched between the a and /3 subunits. 

The very symmetric arrangement with short distances between the pig- 
ments provides a new impulse to the discussion about, e.g., the exciton trans- 
fer regime and the exciton delocalization in LH2. The extent of the exciton 
delocalization, which could be reduced by dynamic and static disorders, has 
been discussed. The principal question reads: Are the states contributing to 
the optical and transport properties of the LH2 ring localized or delocalized? 

In the past knowledge of the energy transfer was mainly derived from steady 
state spectroscopic experiments leading to several absorption bands at dif- 
ferent wavelengths. Pigment molecules B850, that we are dealing with, are 
characterized by the absorption wavelength in the LH2 subunits of the ASs 
of purple bacteria. At room temperature the solvent and the protein environ- 
ment fluctuate with characteristic time scales ranging from femtoseconds to 
nanoseconds. The dynamical aspects of the system are reflected in the line 
shapes of the electronic transitions. To fully characterize the line shape of a 
transition and thereby the dynamics of the system, one needs to know not 
only the fluctuation amplitude A (coupling strength) but also the time scale 
of each process involved. The observed linewidths reflect the combined influ- 
ence of static disorder and exciton coupling to intermolecular, intramolecular, 
and solvent nuclear motions. The simplest approach is to decompose the line 
proflle into homogeneous and inhomogeneous contributions. In more sophisti- 
cated models, each process is deflned with its characteristic time scale as well 
as a coupling strength 
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Both the exciton-phonon couphng and the static disorder lead to local- 
ization of excitons. One tries to avoid the influence either of the dynamical 
disorder by measuring at very low temperatures or of the static disorder by 
performing single-molecule spectroscopy experiments. The absorption, circu- 
lar dichroism, and hole-burning spectra of the LH2 complex are reproduced 
reasonably well by theoretical treatments, that consider Frenkel exciton states 
of the entire complex but invoke moderate disorder in the excitation energies 
of the individual BChls 0. 

In a homogeneous system, in which the Qy transition dipole moments of the 
B850 BChls all lie approximately in the plane of the ring, most of the dipole 
strength of the B850 band would come from a degenerate pair of orthogonally 
polarized transitions at an energy slightly higher than the transition energy 
of the lowest exciton state. However, energetic or structural disorder that 
disrupts the symmetry of the complex will make both the low-energy transition 
and transitions at higher energies weakly allowed. 

Time- dependent experiments (absorption, fluorescence, photon echo, etc.), 
which used (sub-)picosecond light pulses of low energy with high repetition 
rate tunable through the infrared absorption bands of the various BChl pig- 
ments in a variety of bacterial ASs, made it possible to study the long-time as 
well as femtosecond dynamics of the energy transfer and relaxation []T],[T^-|I1|. 
The nature of the rapid relaxation is of interest because it depends on how 
photosynthetic antenna complexes absorb and transfer energy and because it 
may provide an experimental window to protein dynamics on very short time 



scales 15 



The interpretation of time-dependent experiments on femtosecond time 
scale requires a theory which incorporates both dynamic disorder and differ- 
ent kinds of static disorder. Recently, models based on extended exciton states 
with moderate static disorder A ^ J/2 have found acceptance in the interpre- 
tation of the experiments [^^-[1^, |IB|-^ in which it has been shown that the 
elementary dynamics occurs on a time scale of about 100 fs p!5| , pT| -p3[]. For the 
ring of BChls in the LH2 antenna complex Kumble and Hochstrasser ||2^ has 
presented a time-domain analysis of the effects of local energy inhomogeneity 
upon the dynamics of optical excitations. They examined the manifestation 
of disorder scattering in polarized femtosecond spectroscopy and the degree to 
which exciton delocalization is revealed in emission and transient absorption 
anisotropy measurements. The time evolution of the quantum states, prepared 
by impulsive excitation of a statically disordered circular aggregate model for 
LH2 antenna complexes, have been calculated exactly quantum mechanically 
for varying degrees of inhomogeneity (static disorder), but without taking into 
account the interaction with a heat bath (dynamic disorder). For a Gaussian 
distribution of local energies, the dynamics of coherence loss (scattering) has 
been explored as a function of the ratio of the standard deviation A of the dis- 
tribution to the intersite interaction energy J. (Please note that Kumble and 
Hochstrasser used the symbol a for the standard deviation in their paper ||24|| .) 
It has been found, that modest degrees of disorder (A/ J ^ 0.4 — 0.8) are suffi- 
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cient to cause scattering on a sub- 100 fs time scale. The estimate of static dis- 
order emerging from this investigation should be considered as an upper limit. 
By neglecting the influence of phonons, Kumble and Hochstrasser's model does 
not take into account the dynamic scattering effects which contribute to de- 
phasing of the initial wave packet state and which promote thermalization of 
the dephased populations. 

Nagarjan et al. [^] measured the changes in absorption and stimulated 
emission caused by excitation with 35 fs pulses centered at 875 nm. Substan- 
tial relaxation on the time scale of 10-100 fs and an anomalously large initial 
anisotropy of 0.7 in the bleaching and stimulated-emission signal has been 
observed. These results have been interpreted in the simple model of a homo- 
geneous system in which excitations are delocalized over the whole ring. The 
high initial anisotropy was ascribed to coherent excitation of a degenerate pair 
of states with allowed transitions and the relaxations to electronic dephasing 
and equilibration with states at lower energies which have forbidden transi- 



tions. It has been shown in Ref. []T5[ that excitation of membrane-bound LH2 
complexes with low-intensity femtosecond pulses causes changes in absorption 
and stimulated emission that initially depend on the excitation wavelength but 
relax to a quasi-equilibrium with a time constant of 100 ± 20 fs. The signals 
have an apparent initial anisotropy of approximately 0.5 when the complex 
is excited with broadband pulses, and of 0.35-0.4 with narrower pulses. The 
anisotropy decays to 0.1 with a time constant of about 30 fs. 



Using a density- matrix formalism Nagarjan et al. have shown that the 
initial light-induced signals are consistent with coherent excitation of multiple 
exciton levels in an inhomogeneous ensemble of LH2 complexes and that the 
main features of the spectral relaxations and the anisotropy can be explained 
by electronic dephasing and thermal equilibration within the manifold of ex- 
citon levels. For one initial state three different possibilities of the form of the 
exciton density matrix p have been taken into account. Beside the full matrix 
description, which has been represented as a coherent superposition of the al- 
lowed "one exciton" states, the second is obtained by setting the off-diagonal 
terms of p to zero, which gives an incoherent mixture of one-exciton states 
with the same initial populations. Finally, in the third form the diagonal el- 
ements of p are replaced by the populations for a Boltzmann equilibrium at 
temperature T (with the off-diagonal terms zero). 

Nagarjan et al. [0 have discussed several different processes which could 
contribute to the spectral relaxations and the rapid decay of the anisotropy. 
They concluded, that the main features of the spectral relaxations and the 
decay of anisotropy are reproduced well by a model that considers decay pro- 
cesses of electronic coherences within the manifold of the exciton states and 
the thermal equilibration among excitonic states. They have not attempted to 
calculate the dynamics theoretically. But the model described in [15| appears 
capable of rationalizing the observation that the anisotropy decay is faster than 
the spectral relaxation. 

Moderate static disorder, which leads to electronic dephasing, causes a 
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large decrease in the calculated anisotropy, whereas the relaxation of the over- 
all signal occurs mainly during thermalization, because the two most strongly 
allowed exciton states tend to be close together in energy and to have approx- 
imately orthogonal transition dipoles. Nagarjan et al. [jl5| concluded that the 



static inhomogeneity assumed in their LH2 model would cause the off-diagonal 
terms Pab and pf,a of the density matrix to decay with a time constant of the 
order of 160 fs. This suggests that the more rapid dephasing indicated by the 
anisotropy decay (r ^ 30 fs) is driven mainly by dynamic processes. 

The aim of the present paper is to extend the investigation by Kumble 
and Hochstrasser and also by Nagarjan et al. |15] taking into account the 
simultaneous influence of static and dynamic disorders on the exciton scat- 
tering after impulse excitation and on the time-dependent optical anisotropy. 
Different time dependent spectroscopic properties have been calculated before 



but within the secular approximation [25 



The remainder of this article is organized as follows: In the next section 
we present the model we have used. In section 3 several different theories 
giving the dynamical equations for the exciton density matrix are reviewed. 
Our results are presented in section 4 and in the last section we discuss them 
and draw some conclusions. 



2 Model 

2.1 Hamiltonian of the ideal rings 

In this paper we deal with just one exciton, already created at time t = 0, on a 
ring interacting with a heat bath. The Hamiltonian therefore consists of three 
parts: 

m,n q W m q 

= -^ex + Hph + Hcx-ph- (1) 

The first term is the Hamiltonian H^^ belonging to the single exciton in the 
ideal ring. Operators aj„ or create or annihilate an exciton at site m. 
Jmn for m 7^ n is a transfer integral (also called resonance integral) between 
sites m and n. The diagonal elements Jnn are the local energies e„ at site n 
which could take into account, e.g., a dimerization e2n 7^ ^2n+i- The second 
term, Hph-, describes the phonons. We assume independent heat baths for 
each chromophore and the harmonic approximation for the phonons. 6j and 
hq denote phonon creation and annihilation operators, respectively. The last 
term, He^-ph, represents the interaction between the exciton and the bath. 
The exciton-phonon interaction is assumed to be site-diagonal and linear in 
the lattice displacements. The term denotes the exciton-phonon coupling 
constant. 
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In the framework of the extended exciton states it is assumed that both 
intra-dimer J12 and inter-dimer J23 transfer integrals are strong enough to 
build exciton states spanning the whole rings. Inside one ring the pure exciton 
Hamiltonian H^^ (jl|) can be diagonalized using the wave vector representation 
with corresponding delocalized "Bloch" states and energies. Considering, e.g., 
only nearest neighbor transfer matrix elements Jmn = —JiSm,n+i + Sm,n-i), the 
same local energies e„ at every site and using Fourier transformed excitonic 
operators (Bloch representation) 

afc = E«-e^'" ' f' = JJ^ ' / = 0,±l,...±iV/2 (2) 

the simplest exciton Hamiltonian in k representation reads 

H',^ = J2E,ata,. (3) 

k 

The dispersion of the excitonic energies is given by 

Ek = -2J cos k . (4) 

A splitting of the local degenerate exciton energies in one ring, to a band of 
energies Ek, corresponding to the exciton eigenstates, is produced. 



2.2 Hamiltonian of the static disorder 

The quasistatic disorder is described by the time-dependent part of the Hamil- 
tonian H = + Hi{t). This part 

Hi(t) ='^en(t)a^an (5) 

n 

describes fluctuations of the local excitation energies en(t) due to slow mo- 
tion of the protein environment. Fluctuations of the local excitation energies 
en{t) have exponentially decaying correlation functions. The mean values and 
correlation functions are given by 

= 0, 

(e„(t)em(r)) = 5m„A^ exp(-A(t - r)). (6) 

In what follows the simplification is used, that we work only with pure static 
disorder, limA 0, and a Gaussian distribution for the uncorrelated local 
energy fluctuations e,^ with a standard deviation A. 



2.3 Microscopic parameters 



Despite the long history of theoretical and computer modeling of the experi- 
mental results no unified set of microscopic parameters entering the Hamilto- 
nian H has been revealed up to now [p!0Hl^, pi6H20|, p6 . 
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The one-exciton band of the ideal ring E{k) consists of two groups of 9 tran- 
sitions arranged asymmetrically relative to the zero local site energy. These 
Davydov manifolds appear because the 18 BChl molecules in the B850 ring 
are grouped into heterodimers with different coupling energies within the el- 
ementary dimers J12 and between them J23 as well as different local energies 
ei 7^ £2- 

In the presence of static disorder (A 7^ in ifi) virtually all exciton transi- 
tions gain strength at the expense of the low energy k = ±1 transitions, which 
dominate the spectrum of the ideal B850 ring. There is still a discrepancy 



between the measurement by Freiberg et al. [11] and results of hole-burning 
experiments by Small (see references in []T|). In the former experiments con- 
siderable intensity of about 13% (averaged over a large ensemble) attained by 
the lowest k = transition has been revealed. In interpreting these experi- 
ments it was concluded |jll[ that the shape of the linear absorption spectrum 
of B850 at 8 K could be well reproduced by taking A ^ 0.7. The hole burning 
results at 4 K imply much smaller disorder leading to a relative intensity of 
the k = subband of only 3%. The model parameters obtained by Freiberg 
from fitting the simulated linear and nonlinear absorption spectra to the low 
temperature experimental data are the transfer integrals J12 = 375 cm~^ and 
Ji3 = 20 — 30 cm~^ as well as the strength of the static disorder A = 216 cm~^ 
with FWHM = 507 cm~^. The dynamic disorder leads to a broadening of the 
exciton levels E{k) with a FWHM = 58 cm"^ = J/6.5. 



Novoderezhkin et al. |[T^, ^ fitted the pump-probe spectra of the LH2 



antenna. Taking into account the difference between the local BChl energies 
Ae = ei — €2 in the basic dimer which leads to the splitting of the exciton 
energy band to two subbands, the following set of microscopic parameters : 
the transfer integrals J12 = 400 cm~^, J23 = 290 cm~^, the the homogeneous 
line widths of exciton states in the lower subband Fi^^ = 240 cm~^ and in 
the higher subband Tin = 340 cm~^ and the strength of the static disorder 
A = 450 cm~^ has been used. Koolhaas et al. [1^ have recently obtained the 
transfer integrals J12 = 396 cm~^ and J23 = 300 cm~^. 

In what follows we will use the transfer integral J12 as energy unit. We 
choose J23 = 0.7, Ae = ei — €2 = J12, and A between 0.2 and 1. To convert 
the time r in these units into seconds one has to divide r by 2ttcJi2 with c 
the speed of light in cm/s and J12 in units of cm^^. So r = 1 corresponds to 
21.2 fs for J12 = 250 cm'^ and to 13.3 fs for Ju = 400 cm'^ 



2.4 Nature of the initially prepared state 

For a coplanar arrangement of site transition moments /2„ dip ole- allowed tran- 
sitions from the ground state populate only the degenerate k = ±1 levels of 
the ideal (without dynamic and static disorders) ring with Hamiltonian H^^ 
and eigenstates \k). If static disorder of the local energies is present (A 7^ 0), 
the stationary states, i.e. the eigenstates \a) of the Hamiltonian H^^ + Hi, 
correspond to mixtures of \k) and an excitation will prepare a superposition 
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of the \k) states. 

However, Kumble and Hochstrasser concluded, that in the case of 
pump-pulse excitation the dipole strength is simply redistributed among the 
exciton levels due to disorder. So the amplitudes of site excitations and the 
phase relationships in the initial state are necessarily identical to that of an 
equal superposition of A; = ±1 excitons of the ideal ring. Thus, generally, 
the excitation with a pump pulse of sufficiently wide spectral bandwidth will 
always prepare the same initial state, irrespective of the actual eigenstates 
of the system. The nature of this initial state is entirely determined by the 
selection rules of the unperturbed system. 

We shall use the following definitions for the dipole strength 

N 

f^a = Y. <f^n (7) 
n=l 

of state I a) of the real system and for the dipole strength 

N 

f^a = Yl Cn/^n (8) 
n=l 

of state I a) of the ideal crystal. The coefficients and are the expansion 
coefficients of the eigenstates of the ideal and disordered aggregates in site 
representation. 

The pump pulse is characterized by the polarization unit vector Cj. Then 
the initial condition for the density matrix by pulse excitation with the polar- 
ization Cj is given by (Eq. (la) in [p!5|): 



Pap{t = 0]ei) = ^{ei ■ fla){ei ■ flis), (9) 



where A = J2a{(^i ' i^a){(^i ■ j^a)- This initial condition for pulse excitation |2| 
should be expressed in the eigenstate basis of the real system: 

pabit = 0;ei) = ^^(a|a)pc./3(t = 0;e;)(/3|6) 

a p 

= Mi^i- Pp){a\(^){(3\b), (10) 



^From ([lOD one obtains 



y iKi^^-i^.) ^^(el . ;2„)(el ■ /2^)(a|a)(/3|6). (11) 



A 
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2.5 Static disorder-induced scattering of the initially 
pure state 

Pump-probe and spontaneous emission measurements on the femtosecond time 
scale have revealed ultrafast depolarization in the core and peripheral ASs with 
biphasic kinetics: a fast sub-100 fs component and a slower decay on a 200- 
400 fs time scale were observed in all cases illl. Kumble and Hochstrasser 123 



thoroughly investigated the time dependence of the disorder-induced scattering 
of the initial state |fc ± 1). The evolution of the initially delocalized exciton 

P.(t) = |(fc|e-^^*|fc)r (12) 

has been followed considering only the influence of a static uncorrelated dis- 
tribution of the local site energies in a circular aggregate. 

In the basis \k) of the ideal ring the time evolution can be described as 
scattering of the initially prepared exciton \k = ±1) into those levels to which it 
is coupled by the presence of disorder. The nonexponential nature of scattering 
dynamics has been demonstrated: for A = J/3.5 the initial state decays in 
two distinct phases; a dominant ~ 100 fs component and a slower (200-300 fs) 
decay are evident. 

2.6 Consequences of the static disorder-induced scat- 
tering in optical anisotropy measurements 

Optical anisotropy measurements provide a convenient mean to follow experi- 
mentally the decay of a virtual exciton within the LH2 ring. Since this initial 
state is prepared as a superposition of eigenstates of the aggregate, polarized 
pump-probe (or fluorescence) signals contain contributions from interference 
terms which last so long that the phase coherence between different levels is 
maintained. The decay of these cross terms directly reflects the time scale of 
electronic dephasing. In the present case it therefore represents the survival 
time of a virtual exciton within the disordered system. 



Kumble and Hochstrasser [0] calculated the usual time- dependent anisotropy 
using 

/,\ _ {Sxxjt)) — {Sxyjt)) 

{Sxx{t))+2{Sxy{t)) ^ ' 

where the brackets () denote the ensemble average and the average over the 
direction of the laser pulses with fixed relative directions el and 6*2. For ex- 
citation and probe pulses which are polarized along the and Cy directions, 
respectively, the transient gain signal from the delocalized initial state is given 
by H 

{Sxy{t)) = (I E (e. ■ f^aWy ■ /2,)Cde— ) (14) 

where the indices a and / label the eigenstates of the virtual and disordered 
aggregates, respectively. The index n represents the individual site excitations. 
The gain signal {Sxx(t)) is defined accordingly. 
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In Kumble and Hochstrasser's modeling, the decay of the initial state has 
been directly manifested itself as a depolarization reducing the anisotropy from 
0.7 to 0.3 — 0.35 and a subsequent recurrence to establish a final value of 0.4. 
Kumble and Hochstrasser supposed that equilibration between eigenstates of 
the aggregate would lead to a further decrease of the anisotropy to ~ 0.1, but 



this process has not been considered in their paper |2J]. They estimated, by 
comparison of their results with experiment, that the static disorder in a LH2 
ring has a value of A 0.4 — 0.8 which should be considered as an upper 
limit since dynamic variations of the relative eigenstate energies occurring on 
a 100 fs time scale would contribute to the scattering process. 

2.7 Effects of the dynamical disorder 

As pointed out in Sec. 1 the aim of the present paper is to extend the above 
described theory given in Ref. |Q by including the effect of dynamic scattering 



which contributes to dephasing of the initial wave packet and promotes a ther- 
malization of the dephased populations. To deal with the dynamic disorder, 
one works within the exciton density matrix formalism instead of using only 
the exciton wave functions. 

Taken only those terms into account which are directly connected to the 
polarization the time-resolved fluorescence Pif(a;,t) with initial Cj and final 
polarization e} of the light is given (cf. Chapter 3, Eq. (2.58) in |^) in our 
notation as 

Pif{uj,t) = AJ2T,Pii'{i){ef ■ fli'){ef ■ fli)[S{uj - w/zq) + 5{uj - uio)] = 
I I' 

AEEpii'{t)EE{ef ■ fikWf ■ fij){j\l){l'\k)[5{uo - uovo) + 6{uo - uio)] (15) 

IV k j 

and in case of no interaction with a bath as 

Pif{uj,t) = 

El El' E.E/3e-^*(^'-^")(el ■ PaM ■ flp){l\a){mief ■ pvWf ■ jli) x 

[5{u - ujvq) + 5{uj - uio)]. (16) 

This expression can be directly rewritten (after integration over uj) in the form 
of Kumble and Hochstrasser: 

Sij{t) = J Pif{uj,t)duj = 

I EEEEe-^^^^^'-^^'Hei ■ fi.w. ■ ■ iivm ■ ^i) 

•' I I' a p 

X [5(UJ — OJiiq) + 5{UJ — UJiQ)\duJ = 

/EEEEe-*^^-^"^x 

I I' a B 
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[ei-pio){ei-pip)^^{l\n){n\a){l3\n'){n'\l'){ef-^ 

n n' 



X 



a I n 



/3 V 



J2 (el ■ fla){ef ■ /iO(a|n)(n|/)e-^-'f . (17) 



a,Ln 



We see, that the crucial quantity entering the time dependence of the 
anisotropy is the exciton density matrix. In the next section we shall therefore 
review several methods leading to dynamic equations for the exciton density 
matrix which include microscopically the interaction of the exciton with the 
bath. They could be obtained from, e.g., the general Tokuyama-Mori approach 
only by serious approximations. One of them is the Redfield model with and 



without the so-called secular approximation [28, 29], which has been applied to 



exciton transfer and relaxation problems in AS rings several times. It has been 



shown pO| , pT| that in contrast to the justification of the secular approximation, 
the non-secular terms do not average out in time. 

Our result, presented in section 4, have therefore been obtained with the 
Redfield model without secular approximation. 



3 Dynamic equations for the reduced exciton 
density matrix 

The bath-related relaxation can be described in a variety of ways. Among 
others these are the path integral methods [^Q, the semigroup methods 
36[| , and the reduced density matrix (RDM) theory |3^. Here we concentrate 



on the latter. 

Having defined the Hamiltonian H of the system, one exciton and a bath, 
one has to solve the equations of motion for the complete density matrix a, 
namely the Liouville equation 

d 

ih—ait) = La. (18) 
dt ^ ^ ^ ^ 

This is not a simple task because one has to find the time development of 
all matrix elements of the density matrix - diagonal and off-diagonal - in any 
representation which takes into account the exciton and phonon states. 

Nevertheless, this treatment is not necessary in many cases. Information, 
which is used in further investigations, is in many cases limited. For example, 
the site occupation probabilities Pm{t) of the exciton (the diagonal matrix el- 
ements in site representation after averaging over the bath variables) are in 
many cases the most interesting quantities in theoretical investigations of the 
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exciton transfer. The transfer problem of, e.g. a Frenkel exciton interacting 
with a bath, is usually treated by projecting out the bath degrees of free- 
dom. This yields dynamic equations for the single-exciton density matrix p 
complementing the usual single-exciton Liouville equation by terms describing 
the bath influence. The latter terms must, on the other hand, be properly 
parametrized. 



3.1 Stochastic theories 

In the stochastic treatment of the exciton interaction with the dynamic and 
static disorders the influence of disorder is described by a stochastic process, i.e. 
the interaction part of the microscopic Hamiltonian is replaced by a stochas- 



tic time-dependent model Hamiltonian. In the original version |38| of the 
stochastic Liouville equation method (SLE) the bare exciton is influenced by 
a stochastic field 

H2{t) = hmn{t) atiO-n- (19) 



Grover and Silbey ||39[ suggested that one could understand the excitation 
in molecular aggregates also as a fully dressed exciton (polaron) obtained from 
the microscopic Hamiltonian using a canonical transformation. The stochastic 
field acting on the dressed excitation has a different microscopic meaning as 
that one acting on the bare exciton The stochastic field obtained from 
the microscopic Hamiltonian for the local and linear exciton-phonon interac- 
tion leads 0] , in the bare exciton representation, mainly to local energy fluc- 



tuations while the corresponding stochastic field acting on the polaron leads 
mainly to fluctuations of the renormalized transfer integrals. 

The exciton-phonon interaction in the ring subunits LHl and LH2 of the 
purple bacteria should not be very strong Then the Hamiltonian H2{t) 
models the influence of dynamic disorder mainly via fluctuations of the local 
exciton energies hmmif)- Mean values and correlation functions are given by 

{hnn(t)hmm{r)) = cxp (- A |t - T | ) . (20) 

Multitime correlation functions entering the equations of motion can be cal- 
culated from two-time correlation functions ( pO]) according to different rules. 
Mostly dichotomic, Gaussian and white noise statistics have been applied in 
the past. 

The simplest description of the exciton dynamics in the framework of the 
stochastic theories is obtained in the so-called white noise limit A^/A —>■ 70 
(A 00, A 00) in which the time correlation functions become 6 functions. 
The SLE for the bare exciton density matrix in the Haken-Strobl-Reineker 
parametrization (HSR-SLE) (for a review, see [^) have been used very often 
to describe the coupled coherent and incoherent regime of exciton transfer in 
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molecular aggregates using a few phenomenological parameters 7: 



d_ 
di 



2r 



-^([^o,p(t)])^„ 

it)] 

P 

)[2r 



(21) 



In most applications of the HSR-SLE nearest neighbor parameters 70 = 7nm 
for n = m, 71 = 7„m for n = m± 1 and 71 = 7nm for n = m± 1 have been used. 
We have shown PDI that in case of a weak interaction which is site-diagonal 
and linear in the lattice displacements only 70 7^ 0. Such treatment has been 
used by us |^2|-^ and Leegwater |^ in the investigation of exciton transfer 
in ring structures. 



3.2 Redfield theory 

Provided that the excitation dynamics is not very fast, the coupling is rather 
weak, the bath and the system are initially statistically independent, and ex- 
cept for the initial time- interval t ~ (t^ = dephasing time of the bath), the 
approximate equation turns to the form first derived by Redfield p8|,|2l 



ij^p{t) = ^[H,p{t)]+np{t), (22) 

where the Redfield relaxation superoperator TZ (as a linear time-independent 
superoperator) describes the influence of dynamics of the medium and its in- 
teraction with the exciton. In the modern fashion, derivation of Eq. (|22|) is 
perhaps best presented in or In the basis of eigenstates \u)^ . . . 



of the relevant system part the Redfield theory reads 
d i 

'l^Pf^^i^) = -^E^^p^^{t) + T^tivi^'u'Pi^'u'it) (23) 

where E^j^u = Efj_ — = hu^i,. The physical interpretation of the so-called 
Redfield tensor IZ^^^iyi, is straightforward: Each element is essentially a rate 
constant that relates the rate of change of reduced density matrix elements 
p^iy to the current value of the other matrix element p^'i.'. Thus 7^^/^/^/^ is the 
relaxation rate of the depopulation of eigenstate while TZf^/^uu represents the 
rate of population transfer from state |z/) to the state The elements TZ^ufiu 
represent the dephasing rates of the off-diagonal density matrix elements p^^. 
An important property of the Redfield tensor is that when the bath is treated 
quantum mechanically, the physically required detailed balance conditions are 
naturally satisfied, so that the subsystem relaxes to thermal equilibrium with 
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the bath. This condition means that the ratio of forward and reverse rate 
constants between states and |z/) is 

1^ = eM-pnu;,.) ■ (24) 

Other elements of the Redfield tensor cannot be given such a simple in- 
terpretation. Keeping them as just fitting semi-empirical parameters leads to 
technical complications as the total number of elements of TZ increases as A^^ 
(A^ is the number of eigenstates involved, i.e. two for the dimer). Hence, the 
idea is to omit those elements of TZ which are physically less relevant provided 
that this omission does not appreciably influence the solution. That is why al- 
ready Redfield himself p8|,|^ suggested the so-called secular approximation. 
This approximation amounts to neglecting those elements of TZ that couple 
elements of the density matrix Pf^i^/ with incommensurate frequencies E^y/h 
determined in the first term on the right hand side of (E^). In other words, it 



means to put R^ufi'u' = whenever E^^ ^ E^iyi. It has been shown pO| , pl| that 
this approximation leads to non-physical effects in the transfer processes and is 
thus, in spite of its apparent success in many cases, physically not acceptable if 
we discuss, e.g., relaxation among states other than eigenstates of the Hamil- 
tonian of the system. Other approximations to the Redfield theory, which have 
not been applied here, have been used and tested in the context of electron 
transfer [^,^. Aspects concerning the efficient numerical implementation of 



the Redfield equation are discussed elsewhere . 

Although being aware that Redfield theory is limited to the high temper- 
ature regime i.e. the system relaxation times have to be large compared to 
the thermal time h/{ki,T) we assume here here that this theory gives reason- 
able results also at low temperatures. In the present study Redfield theory is 
applied to the full temperature range. 

3.3 Generalizations of the HSR-SLE method in the local 
representation 



The original HSR-SLE method |Q represents the simplest stochastic treat- 
ment of the microscopic interaction Hamiltonian, but the application of the 
white noise is very restrictive. It leads to improper long time asymptotics of 
the occupation probabilities of the extended exciton states Pfc(oo). Many at- 
tempts have been made to determine more general equations using different 
representations of the exciton interacting with the bath and different ways of 
obtaining the dynamic equations for the exciton density matrix from the full 
Hamiltonian (|l]). The equations obtained in such a way could, after some fur- 
ther approximations and manipulations, be rewritten |^8|-pO[] in a form which 



resembles the nonconvolutional SLE but brings several important improve- 
ments. 

Capek P8|-|50| applied several different ways of obtaining of the convolu- 



tional and convolutionless dynamical equations for the exciton density matrix. 
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The equations of motion obtained read 

d 
di 



^Pmn{t) = YllUJmn,pq{t)ppqit), (25) 
PI 



where 

When the typical system times arc much longer than the bath correlation time, 
i.e. in the Markovian limit, one obtain a system of equations with constant 
coefficients 

^mn,pq '^\,^mp^nq{^n ^m) ~l~ Jqn^mp Jmp^n<^i 

jXQ ^ -S i« -S * 

''^^''mn,pq — ^mp'^nm ^nq'^mn ' 



^'mn = '^T.^l{G^,-G%)Y.GlY.{^2\r){r\v,){v,\m){p\v,)x 

k r 1/1,1/2 

\ E^^ - Ey^ + hwk + is E^^ — E^2 - hujk + is 

Here |z/) and Ey label the corresponding eigenvalues and eigenvectors of H^^. 
Capek suggested the parametrization 

Al,^ = Re = Im (28) 

In the fast carrier regime the coefficients become negligible (no renormal- 
ization of Jmn) and can be omitted. 

The play the same role as the coefficient 70 = 7mm of the HSR 
parametrization. The coefficients provide the bath-induced coupling of 
to the off-diagonal to diagonal elements of Pmn(^)- 




4 Results 

4.1 Equivalence of Redfield and Capek's equations 

We have used the dynamical equations of motion for the reduced cxciton den- 
sity matrix in the Capek's form (25) and (26) to express the time dependence 
of the optical properties of the model LH2 ring of BChls in the fs time range. 

For linear local exciton-phonon coupling the Hamiltonian H^x-ph can be 
written as 

Hex—ph ^ ] ^m^m-i (^^) 
m 

where the system operators K^^ and bath operators $m are given by 

% 
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If we express the equation of motion for the reduced density matrix in 
Born approximation (second order in exciton-bath couphng) in terms of Km 
and $m, and apply the Markov approximation (in the interaction picture), the 
dissipation term looks like 



'dp\ 

^ / diss 



E 



h " JO 

m,n 



(30) 



where Cmn{t) is the bath correlation function 

C„^nit) = ($™(t)<l'„(0))B. (31) 

Here we assume that the expectation values ($m) vanish. The original Redfield 



equations are derived in eigenstate representation [^,|2^. Here we will express 
Eq. (|30|) in site representation 

/ dprs \ 



dt 



diss 



~ E r Cmnir) E {{r\Km\t) mi'\-T)\j)p,,{ty 

m,n''^ I ij 

{r\Ki'\'r)\^){J\KMP^At)] 
-CUr)Y.{{^\Kn^\^)mi'\-r)\s)p,,it)- 

m^'\-r)\j){j\Km\s)pr.{t))}, (32) 



where the indices i, j, r, s denote sites. 
For Km as given above we have 



(33) 



Using the simple expression for the time evolution operator U{t) in energy 
representation (a|f/(r)|/3) = ^q./jc"^'^"'^ one gets 

{r\K^^\~T)\s) = {r\U{r)KmUHr)\s) = 

i,j a,l3 a,l3 

(34) 

In this treatment the matrix elements of K^^ are evaluated exactly. In the 
context of electron transfer the so-called diabatic damping approximation has 
been used and tested PT],ES|. In this approximation the influence of the site 
coupling on dissipation is neglected. 
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After substituting Eqs. ( ^31) and (|34D into (|3^ ) the final form of the equation 
of motion for the reduced density matrix in local site representation reads 



' dpr 



POO 



\ J diss ^ n,j a,l3 

X m.n{r) - CUr)){r\a){a\n){n\[3){mp,s{t) 
-{C:^{r) - C:Sr)) {3\a){a\n){n\P){P\s)pr,{t)} . (35) 

If one now introduces the parameters 
1 r°° 

^Ts = T^T. rfre— '^''-$:(C:„(r)-C;jr))(m|a)(a|n)(n|/?)(/?|r) (36) 

it is possible to express Eq. ( |35D in the same form as it was done by Capek ^0 

( dprs \ 
V dt ),^^^ 

where 

iS^rs^pq ^TpAgj. SsqA^g . (^^) 

So Redfield theory without secular approximation and Capek's theory after 
Markov approximation are equivalent. 

In what follows we have used a simple model for Cmn- Each site (chro- 
mophore) has its own bath with one bath coordinate, and the various coordi- 
nates are completely uncorrelated. Further it is assumed that these baths have 
identical properties Then only one correlation function C(a;) is needed 

Cmn{uj) = 6mnC{Lj), (39) 

with 

Cmnioo) = / C„„(r)e^-^dr = 2Re / C^„(r)e^'^^cir. (40) 

J-oo Jo 

In case of real eigenvectors (which always can be chosen for this system) 
{m\a) and if we omit the imaginary parts of A^g which cause only renormal- 



ization of transfer integrals J |]50|, we can write 
1 
2E 



AT, = ReATs = ^Y.CM{m\a){^\r){{a\s){s\^) - («|r)(r|/5)) (41) 



a/3 

and after introduction of the spectral density J{uj) |^ 

C{uj) = 27Th'^[l + nB{uj)][J{uj) - J{-uj)], (42) 
it is possible to write 

= E [1 + (oOafs)] [J{uJap) " ^(^/3a)] X 

a/3 

{m\a){P\r){{a\s){s\P) - {a\r){r\P)) 
= -n Yl JM{{a\s){s\/3)-{a\r){r\P))x 

af3{e0>£a) 

{[1 + nB(cu,;3)](m|a)(/?|r) + nB{uJc,p){m\(3){a\r)} . (43) 
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In the following we have modeled the spectral density by the form usually 
used dealing with the BCHls in the LH2 ring unit p|, p5| , |45[1 

,2 



J{uj) = &{uj)jo 



2^* 



(44) 



It has its maximum at 2co'c- We shall use (in agreement with [^) jo = 0.2 or 
jo = 0.4 and cUc = 0.2. 



4.2 Disorder-induced scattering of the initial state 

In the unperturbed ring the time evolution of a exciton in the presence of 
static site energy and dynamic disorder leads to a loss of its coherent char- 
acteristics, i.e., destruction of the phase relationships between the individual 
eigenstates comprising the initial wave packet state. We have calculated the 
time dependence of the occupation probability of the initial state after impul- 
sive excitation, namely 

Pk{t) = pkkit) (45) 

for an equal superposition of A; = ±1 states. In the basis of the unperturbed 
system the time evolution can be described as scattering of the initially pre- 
pared exciton into the levels to which it is coupled by the presence of disorder. 
The scattering dynamics of a single ring after averaging over more than 1000 
realizations of the static disorder A is shown in Fig. 1 for four static disorders. 
The influence of the dynamic disorder is displayed for different temperatures 
for three values of dissipation. The curves (i) with no dissipation do not 
coincide with the corresponding results by Kumble and Hochstrasser, Fig. 2 
in Their figures correspond to the time dependence of P+iif) or P_i(t). 
In Fig. 1 one can see a pronounced effect of the dynamic disorder (curves (ii) 
- (v)) on the disorder-induced scattering of the initial state, while the results 
without dynamic disorder (the Redfield tensor set to zero) (i) coincide with 
those of Kumble and Hochstrasser (when comparing the time dependence 
of -P+i(t) or P-i{t) with Fig. 2 in Ref. 0). With dissipation the populations 
equilibrate while without dissipation the populations first decay rapidly but 
then increase again. Without static and dynamic disorder this would lead to 
full revivals. For short times some reminiscences of these revivals are seen if 
one neglects dissipation. For low temperatures (cases (ii), (iv)), the dynamic 
disorder leads to a further decay of the occupation probabilities Pi _i(t) after 
an initial fast decay. Low temperature correspond to ksT/J^ = 0.001 i.e. for 
Ji2 = 250 cm~^ the temperature is T = 0.367^. For room temperature the 
inclusion of dynamic disorder leads to faster initial decay. 



4.3 Consequences of the disorder-induced scattering in 
optical anisotropy measurements 

The time-dependent anisotropy (|13]) has been calculated for excitation pulses 
polarized along the directions (cj = Cx ) and probe pulses which are polarized 
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along the and Cy directions (e/ = Cx and e/ = Cy). It is shown in Fig. 2. 
Similar conclusions as for Fig. 1 can be drawn from Fig. 2. It was concluded in 
Ref. based on measurements by Chachisvilis et al. |^ that the time decay 
of the anisotropy during the first dozens of fs is temperature independent in the 
case of LH2 subunits. Our calculation show that such result can be obtained 



only for A > 0.8. But because the time resolution of the experiments in |5T 



was not too high this very restrictive statement about the strength of the static 
disorder can only be made with caution. We expect that some temperature 
dependence can be seen using experiments with shorter laser pulses. 



5 Discussion and concluding remarks 

We have used the LH2 ring structure of real antenna subunits of purple bac- 
teria photosynthetic systems as a motivation for our investigation of the time 
dependent optical spectra of our molecular aggregate ring model. In a series of 
papers the steady state optical absorption line shape of the Frenkel excitons on 
cyclic model oligomers have been calculated in presence of both dynamic 
and quasistatic disorders. Correlation with Freiberg's experimental data gives 
an estimate for the strength of the static disorder A = 0.7 [ p^ . 

In the present investigation we study the time-dependent optical properties 
under the influence of two simultaneous site energy fluctuations - quasistatic 
and dynamic within the LH2 ring antenna subunits of purple bacteria PSUs. 
Despite a long history of experimental and theoretical investigations no final 
conclusion has been done on the ratio A/J12. Kumble and Hochstrasser con- 
cluded to reach a decay in time-dependent optical experiments below 100 fs 
one needs A k, 0.4 — O.8J12. They did not include dynamic disorder in their 
calculations. We have extended those calculations by including dissipation. 
We were able to show that in this case smaller values of A/ J12 are necessary 
than predicted by Kumble and Hochstrasser. 

Our investigation has been based on several simplifications: 

A) The time development of the reduced exciton density matrix has been 
described by the long time limit of Capek's dynamical equations of motion 
(Markovian limit). As shown in this paper Capek's equations are equivalent 
to the Redfield equation without secular approximation. More precise descrip- 
tion of the time development of the exciton density matrix during first dozens 
of fs requires certainly more advanced description with time dependent coeffi- 
cients in Eq. (25) (non-Markovian treatment). We had discussed such results 
obtained by us only for a dimer using different methods for description of the 
time development of the exciton density matrix [Q. In addition we have as- 
sumed that Redfield theory could give reasonable results for low temperature. 

B) We have used the ring symmetry to describe the polarization proper- 
ties of the Frenkel exciton states. Recent results obtained by single molecule 
spectroscopy can only be interpreted admitting the presence of a C2 



distortion of the ring. It has, up to now, not been concluded whether such 
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6*2 distortion of the LH2 ring is present also in samples in vivo. Calculations 
which take into account the possible C2 distortion of the ring are on the way. 

C) Kumble and Hochstrasser have pointed out that the information needed 
for a complete description of the exciton bath interaction is not available for 
the LH2 units. We have worked with a form of the spectral density which has 
been commonly used previously Work on a better description of the 



exciton-phonon interaction model is in progress. 
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Figure captions: 



Fig. 1: Time dependence of the occupation probability (P±i) for static disorder 
A = 0.4 (a), A = 0.6 (b), A = 0.8 (c) and A = 1.0 (d). The influence of the 
dynamic disorder is displayed by curves for low (ii), (iv) and room (iii), (v) 
temperature for jo = 0.2 (ii), (iii), and jo = 0.4 (iv), (v) compared to jo = 0.0 
(i)- 



25 




26 




27 




28 




29 



Fig. 2: Same as in Fig. 1 but for the anisotropy r{t). 
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